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Abstract — The  Nyquist  stability  criterion  is  extended  to  a 
class  of  spatially  periodic  systems  with  spatially  distributed 
inputs  and  outputs.  It  is  demonstrated  that  the  exponential 
stability  of  this  class  of  systems  can  be  guaranteed  by  checking 
the  Nyquist  stability  criterion  for  a  family  of  finite-dimensional 
systems.  In  order  to  show  this  result,  a  new  version  of  the 
argument  principle  is  derived  that  is  applicable  to  systems  with 
infinite-dimensional  input/output  spaces  and  unbounded  system 
operators. 

1.  Introduction 

Of  the  frequency-domain  methods  of  stability  analysis,  the 
Nyquist  criterion  is  of  particular  interest  as  it  offers  a  simple 
visual  test  to  determine  the  stability  of  a  closed-loop  system 
for  a  family  of  feedback  gains  [1]. 

Extensions  of  the  Nyquist  stability  criterion  exist  for 
certain  classes  of  time  periodic  [2]  and  distributed  systems 
[1],  where  by  distributed  it  is  meant  that  the  state-space 
of  the  system  is  infinite-dimensional.  Time-delay  (retarded) 
systems  and  systems  governed  by  partial  differential  equa¬ 
tions  (spatially  extended  systems)  are  examples  of  distributed 
systems.  But  to  the  authors’  best  knowledge,  the  existing 
literature  on  the  Nyquist  stability  criterion  for  distributed 
systems  deals  only  with  those  systms  that  have  finite¬ 
dimensional  input  and  output  spaces.  An  example  of  such 
a  system  would  be  one  described  by  a  partial  differential 
equation  with  pointwise  sensing  and  pointwise  actuation. 

In  contrast,  in  the  present  work  we  aim  to  extend  the 
Nyquist  stability  criterion  to  a  class  of  spatially  periodic  sys¬ 
tems  that  possess  spatio-temporal  (i.e.  infinite-dimensional) 
input  and  output  spaces,  and  possibly  unbounded  system 
operators. 

Some  of  the  difficulties  in  the  application  of  the  Nyquist 
criterion  to  the  class  of  systems  described  above  include  (a) 
the  concept  of  the  characteristic  function  and  its  zeros  are 
not  immediately  extendable  to  unbounded  system  operators, 
and  (b)  due  to  the  infinite  dimensionality  of  the  input/output 
spaces  one  has  to  locate  an  infinite  number  of  eigenloci. 

Hence  the  main  contributions  of  this  paper  are  to  (a)  derive 
a  version  of  the  argument  principle  that  lends  itself  to  un¬ 
bounded  system  operators,  and  (b)  show  that  one  can  always 
truncate  the  infinite-dimensional  system  operators  so  that  the 
problem  can  be  reduced  to  checking  the  Nyquist  stability 
criterion  for  (a  family  of)  finite-dimensional  systems. 
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Our  presentation  is  organized  as  follows:  Section  II  gives 
a  short  overview  of  the  frequency-domain  representation 
of  periodic  operators.  We  lay  out  the  problem  setup  in 
Section  III  and  describe  the  general  conditions  for  stability  of 
spatially-periodic  systems  in  Section  IV.  Section  V  contains 
the  main  contributions  of  the  paper,  where  the  Nyquist 
stability  criterion  is  developed  for  spatially-periodic  systems. 
We  apply  the  theory  to  a  simple  example  in  Section  VI, 
and  finish  with  some  conclusions  and  directions  for  future 
research  in  Section  VII.  All  proofs  and  technical  details  have 
been  placed  in  the  Appendix. 

Notation:  A:  €  M  denotes  the  spatial-frequency  variable, 
also  known  as  the  wave-number.  S(r)  is  the  spectrum  of 
T,  Sp(r)  its  point  spectrum,  and  p{T)  its  resolvent  set. 
(t„(T)  is  the  n*  singular- value  of  T.  denotes  the 

bounded  operators  on  the  compact  operators  on 

and  the  nuclear  operators  on  i.e.  operators 

T  that  have  the  property  <  oo;  C 

C  ££§{£^).  tr[T]  denotes  the  trace  of  T  and  det[T] 
its  determinant.  C+  and  C”  denote  the  closed  right-half  and 
the  open  left-half  of  the  complex  plane,  respectively,  and 
j  :=  v/~l-  C{zo',  ‘P)  is  the  number  of  counter-clockwise 
encirclements  of  the  point  2:0  G  C  by  the  closed  path  ip. 

II.  Preliminaries 

Let  us  consider  a  linear  time  invariant  (LTI)  spatially 
distributed  system  of  the  form  dttp  =  Afi,  where  A  is  a 
spatial  operator  on  L^(IR)  and  is  a  spatio-temporal 

function.  An  example  of  such  a  system  would  be  the  heat 
equation  on  the  real  line  dtfi  =  a;  G  M. 

Now  for  those  systems  where  A  is  a  spatially  invariant 
operator  (i.e.,  the  action  of  A  on  i/;  can  be  represented  by  a 
convolution),  Fourier  methods  can  significantly  simplify  the 
problem.  For  example  in  the  heat  equation  introduced  above, 
is  a  spatially  invariant  operator  on  L^(IR),  and  thus  the 
system  can  be  rewritten  in  the  spatial-frequency  domain  as 
dtfi  =  Here  k)  is  the  spatial  Fourier  transform 

of  A{k)  =  {jkY  =  —k'^  is  the  Fourier  symbol  of 

the  spatial  operator  d^,  and  fc  G  M  is  the  spatial-frequency 
variable. 

In  general  if  an  LTI  system  is  composed  of  only  spatially 
invariant  operators  then  as  illustrated  above,  the  spatial 
Fourier  transform  ‘diagonalizes’  all  such  operators  (i.e.,  turns 
them  into  multiplication  by  a  function  of  k).  Therefore  the 
Fourier-transformed  system  collapses  to  a  ‘continuum’  of 
finite-dimensional  LTI  systems  parameterized  by  k.  This 
means  that  the  original  infinite-dimensional  problem  has 
been  effectively  ‘decoupled’  in  the  spatial-frequency  domain. 
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But  in  this  paper  we  will  be  dealing  with  the  broader  class 
of  spatially  periodic  systems,  i.e.,  systems  of  the  form  dtip  = 
A-ip  where  A  is  a  spatially  periodic  operator  on  L^(IR).  An 
example  of  such  a  system  would  be  dt'tp  =  +  cos{x)'tp. 

In  light  of  the  previous  discussion,  the  difficulty  faced 
in  the  analysis  of  spatially  periodic  systems  becomes  clear. 
For  the  Fourier  methods  described  above  to  simplify  the 
representation  of  a  spatially  periodic  system,  the  Fourier 
transform  has  to  simultaneously  diagonalize  the  spatially  in¬ 
variant  and  spatially  periodic  operators  present  in  the  system. 
But  this  is  not  possible  since  these  two  classes  of  operators, 
in  general,  do  not  commute.'  For  example  take  the  two 
spatial  operators  dx  and  cos(x).  Clearly  cos(a;)  dx  P’ix)  ^ 
dx  cos(a;)  tp{x)  =  —  sin(a:)  ipix)  +  cos(x)  dx  tp{x),  for  a 
general  nonzero  differentiable  function  ip{x).  Thus  cos(x), 
viewed  as  a  spatially  periodic  pure  multiplication  operator, 
does  not  commute  with  the  spatially  invariant  differentiation 
operator  dx- 

Yet  this  does  not  mean  that  Fourier  methods  are  not 
applicable  or  are  not  useful  in  the  analysis  of  periodic 
systems.  Quite  to  the  contrary,  analysis  of  periodic  systems 
in  the  Fourier  (frequency)  domain  often  leads  to  numerical 
tractability  and  valuable  insight. 

In  [3]  and  [4]  it  is  demonstrated  that  a  general  spatially 
periodic  operator  A  on  L^(IR)  with  spatial-period  X  =  27r/n 
can  be  represented  in  the  Fourier  domain  by  a  family  of 
operators  on  These  operators  can  be  written  as  bi-infinite 
matrices  Ae  parametrized  by  a  variable  6  G  [0,12].  In  this 
paper  we  only  give  the  Ae  representation  of  certain  special 
subclasses  of  spatially  periodic  operators,  namely  spatially 
invariant  operators  and  spatially  periodic  pure  multiplication 
operators,  with  the  additional  note  that  any  spatially  periodic 
operator  can  be  written  as  the  summation  and/or  multiplica¬ 
tion  of  a  countable  number  of  such  operators. 

1)  A  spatially-invariant  operator  A  has  the  diagonal 
representation 


Ae  = 


A{9  -\-  fin) 


where  A(  • )  is  the  Fourier  symbol  of  A,  and  for  every 
given  0  G  [0, 12]  the  diagonal  elements  of  Ae  are  the 
equally-spaced  samples  {A{0  +  nl2)}„gz  of  A{  ■ ). 

2)  A  spatially-periodic  pure  multiplication  operator  F  has 
the  Toeplitz  representation 


So 

/-I 

f-2 

/i 

So 

f-1 

/2 

fl 

So 

’TS  =  ST  is  a  necessary  and  sufficient  condition  for  T  and  S  to  be 
simultaneously  diagonalizable  by  a  transformation  U. 


where  //  are  the  Fourier  series  coefficients  of  F{x), 
i.e.  F{x)  =  /t  Notice  that  T  is  inde¬ 

pendent  of  d. 

Remark  1:  If  A  is  a  bounded  operator  on  L^(IR),  then 
suPfcgK  |A(/c)j  <  M  for  some  M  >  0,  and  thus  Ae  is  a 
bounded  operator  on  for  every  0  G  [0,12].  ■ 

III.  Problem  Setup 

Consider  the  spatially-invariant  system  S°* 

{dttp){t,x)  =  {Ail)){t,x)  +  {Bu){t,x), 

y{t,x)  =  lcip){t,x),  (1) 

where  /  G  [0,  oo),  a:  G  M,  A,  B,  and  C,  which  we  call  the 
system  operators,  are  spatially  invariant  and  are  all  defined 
on  a  dense  domain  ^  C  L^(M).  u,  y  and  ip  are  the  spatio- 
temporal  input,  output  and  state  of  the  system,  respectively. 
Clearly,  for  any  given  time  t,  ip{t,  ■)  is  a  spatial  function  on 
L^(M),  and  thus  (1)  is  an  infinite-dimensional  linear  system. 

Next,  we  place  the  spatially-invariant  system  S°'  in  feed¬ 
back  with  a  spatially-periodic  operator  ''yF{x),  ||F'(a;)||  =  1, 
7  G  C,  to  form  the  closed-loop  system  8“^’  as  in  Figure  1 .  It 
is  our  aim  here  to  determine  the  stability  properties  of 
as  the  feedback  gain  7  varies  in  C. 


Fig.  1.  The  closed-loop  system  S' 


For  simplicity,  in  this  paper  we  assume  that  all  operators 
A,  B,  C  and  F  are  scalar.  We  also  make  the  following 
assumptions  on  S°'. 

Assumption  (i):  A  is  such  that  Re(A(/c))  <  (3  for  |fc|  > 
Ki,  and  |A(fc)|  >  for  |fc|  >  K2,  for  some 

/3  <  0,  a  >  0,  77  >  0,  ATi  >  0  and  K2  >  0, 

Assumption  (ii):  B  and  C  are  bounded  operators. 

[3]  and  [4]  show  how  (1)  can  also  be  represented  as 

{dtipe)(t)  =  {Aeipe)it)  +  {Beue){f), 

ye{t)  =  {Ceipe){t),  (2) 

with  6  G  [0,  n],  with  Ae,  Be,  and  Ce  being  the  bi-infinite  ma¬ 
trices  introduced  previously.  (2)  describes  the  open-loop  sys¬ 
tem  Sf  with  temporal  impulse  response  Ge{t)  '■=  Ce  Be, 
and  transfer  function 


Ge{s)  ■=  Ce  {sX  —  Ae)  ^  Be 


C{9  +  Q.n)B{9  +  Cln) 
s  —  A{9  +  fin) 


Finally,  using  the  same  bi-infinite  representation  for  the  pe¬ 
riodic  operator  yF{x)  to  get  yF,  and  placing  it  in  feedback 
with  Sg*,  we  obtain  the  close-loop  system  S^*  in  Figure  2 
with  A-operator  Af  :=  Ae  —  Be  yT Ce- 


Fig.  2.  The  closed-loop  system  Sg. 


IV.  Stability  of  Linear  Spatially-Periodic 
Systems 

A  semigroup  on  a  Hilbert  space  is  called  exponentially 
stable  if  there  exist  constants  M  >  1  and  a  >  0  such 
that  ||e"^*||  <  for  f  >  0.  It  is  well-known  [5]  [6] 

that  if  A  is  an  ifinite-dimensional  operator,  then  in  general 
S(A)  C  C“  is  not  sufficient  for  the  exponential  decay  of 
||e^‘||.  In  this  paper  we  focus  on  systems  which  do  satisfy 
the  so-called  spectrum-determined  growth  condition,  namely 
systems  for  which  S(A)  C  C“  implies  exponential  decay  of 
the  semigroup.  Examples  of  such  semigroups  are  numerous 
and  include  analytic  semigroups  [7]  [8]. 

In  [3]  it  is  shown  that  for  a  general  spatially  periodic 
operator  A  we  have 

S(A)  =  U  S(^e).  (4) 

ee[o,n] 

Thus  to  prove  C  C“,  as  needed  to  guarantee  the 

exponential  stability  of  8“^',  it  is  necessary  and  sufficient  to 
show  that  Ag  =  Ae  —  Bq  Cg  has  spectrum  inside  C“ 
for  all  6  e  [0,  H].  In  the  next  section  we  aim  to  develop  a 
graphical  method  of  checking  whether  or  not  C  C“. 

Also,  henceforth  in  this  paper  wherever  we  use  the  term 
stability  we  mean  exponential  stability. 

V.  The  Nyquist  Stability  Criterion  for  Spatially 
Periodic  Systems 

A.  The  Determinant  Method 


To  motivate  the  development  in  this  section,  let  us  first 
consider  a  finite-dimensional  (multi-input  multi-output)  LTI 
system  G'(s)  placed  in  feedback  with  a  constant  gain  7/. 
In  analyzing  the  closed-loop  stability  of  such  a  system,  we 
are  concerned  with  the  eigenvalues  in  C"*"  of  the  closed-loop 
A-matrix  A‘^’.  If  s  is  an  eigenvalue  of  then  it  satisfies 
det[s/  —  A‘^']  =  0.  Now  to  check  whether  the  equation 
det[s/  —  A‘^']  =  0  has  solutions  inside  C’*',  one  can  apply 
the  argument  principle  to  det[/+7G(s)]  as  s  traverses  some 
curve  X*  enclosing  C^.  More  precisely,  since 


det[/  +  7G(s)] 


det  [si  —  A‘^'] 
det[s/  —  A] 


if  one  knows  the  number  of  unstable  open-loop  poles,  then 
one  can  determine  the  number  of  unstable  closed-loop  poles 
by  looking  at  the  plot  of  det[/  +  7G*('S)] 

But  in  the  case  of  spatially-distributed  systems  the  open- 
loop  and  closed-loop  A-operators  Ag  and  Ag  are,  in  general, 
unbounded.  Hence  it  does  not  make  sense  to  talk  about  the 
characteristic  functions  det  [sX  —  Ag]  and  det  [sX  —  A^],  and 
one  has  to  resort  to  operator  theoretic  arguments  to  relate 


Fig.  3.  The  closed  contour  X)®  traversed  in  the  clockwise  direction  taken 
as  the  Nyquist  path  as  r  — >  oo.  The  indentations  are  arbitrarily  made  to 
avoid  open-loop  modes  on  the  imaginary  axis. 

the  plot  of  det[X  +  to  the  unstable  modes  of 

the  open-loop  and  closed-loop  systems.  But  first  it  has  to  be 
clarified  what  is  meant  by  det\I-\-^TQg{s)]  for  the  infinite¬ 
dimensional  operator  X  +  ^TQg{s).  We  need  the  following 
lemma. 

Lemma  1:  TQe{s)  G  for  all  s  G  p{Ag). 

Proof:  See  Appendix.  ■ 

Since  TQg{s)  G  one  can  now  define  [9]  [10] 

00 

det[X+7X'00(s)]  :=  (l  +  7A®(s)), 

n——oo 

where  A^(s),  n  G  Z,  are  the  eigenvalues  of  Qe{s).  We 
are  now  ready  to  state  a  generalized  form  of  the  argument 
principle  applicable  to  systems  with  unbounded  A-operators. 
Theorem  2:  If  det[X  +  '-^TQq{s)]  f  0  for  all  s  G 

G(0;det[X  +  7X-ge(s)]|^^^g)  = 

tr  /  {sI  —  A'g)~^ds  —  tr  -^  [  {sJ  —  Ag)~^ds  = 

—  (number  of  eigenvalues  of  Ag  in  C^)  + 

(number  of  eigenvalues  of  Ag  in  C^), 

where  D®  is  the  Nyquist  path  shown  in  Figure  3  that  does 
not  pass  through  any  eigenvalues  of  Ag,  and  encloses  a  finite 
number  of  them. 

Proof:  See  Appendix.  ■ 

Remark  2:  V  —  —5^  J^g(sJ  —  Ag)~^ds  is  the  group- 
projection  [11]  [10]  corresponding  to  the  eigenvalues  of  Ag 
inside  X)®,  and  tr[X]  gives  the  total  number  of  such  eigen¬ 
values  [12].  Similarly  tr[— ^  (sJ  —  Ag)~^dsJ  gives  the 
total  number  of  eigenvalues  of  A^*  in  X)^.  Thus  Theorem  2 
allows  us  to  determine  the  number  of  C’*'  eigenvalues  of  Af 
from  knowledge  of  the  number  of  C'*'  eigenvalues  of  Ag 
and  the  number  of  encirclements  of  the  origin  by  the  plot 
of  det[X  +  ^!FQg{s)\  as  s  traverses  X)®.  Notice  that  since 
Ag  =  diag{  •  •  •  ,  A(0  +  fin),  •  •  •  },  the  eigenvalues  of  Ag 
are  known  and  we  have  Sp(Ae)  =  {A{9  +  Lln),n  G  Z}.  ■ 
Remark  3:  It  now  becomes  clear  why  we  have  used  the 
Ag  representation  of  the  operator  A;  in  this  representation 
and  under  the  assumptions  of  Section  III,  the  open-  and 
closed-loop  A-operators  enjoy  the  property  of  having  discrete 


(pure  point)  spectrum.  This  allows  for  an  extension  of  the 
argument  principle  to  be  invoked.  ■ 

Remark  4:  Since  Ae  has  discrete  spectrum,  it  has  no  finite 
accumulation  points  in  the  complex  plane  [12].  In  particular, 
the  eigenvalues  of  Ag  can  not  converge  to  any  finite  point 
jojo  of  the  imaginary  axis.  Yet  this  does  not  rule  out  the 
possibility  of  the  eigenvalues  accumulating  at  ±joo.  But 
Assumption  (i)  guards  against  this  by  requiring  that  A{k)  be 
bounded  away  from  the  imaginary  axis  as  |fc|  ^  oo.  Thus 
the  Nyquist  path  can  be  taken  to  run  to  infinity  along  the 
imaginary  axis  without  any  technical  difficulties.  ■ 

As  a  direct  consequence  of  Theorem  2  we  have 

Theorem  3:  Assume  denotes  the  number  of  eigenval¬ 
ues  of  Ag  inside  C+.  For  S®  taken  as  above,  )  C  C“ 

iff 

(a)  det[X  +  -iTgg{s)]  A  0,  V  s  G  S®, 

and 

(b)  c{()-det[l  +  iTgg{s)]\^^^,)=p%. 

Finally,  the  closed-loop  system  G‘^*  is  exponentially  stable  iff 
C  C-  for  all  e  G  [0,  ft].  ■ 

B.  The  Eigenloci  Method 

The  setback  with  the  method  described  in  the  previous 
paragraph  is  that  to  show  Sp(^g)  C  C“,  Af  =  Ag  — 
Bg^TCg,  for  different  values  of  7,  one  has  to  plot  det[X  + 
for  each  7.  This  also  includes  having  to 
calculate  the  determinant  of  an  infinite  dimensional  matrix. 
This  motivates  the  following  eigenloci  approach  to  Nyquist 
stability  analysis,  which  is  very  similar  to  that  performed  in 
[2]  for  the  case  of  time-periodic  systems. 

Let  A®(s),  n  G  Z,  constitute  the  eigenvalues  of  !Fgg{s). 
Then 

00 

Z  det[I  +  ^Tge{s)]  =  /  (l  +  7A^(s)).  (5) 

n—  —  oo 

But  recall  from  Lemma  1  that  Tgg{s)  G  for  every 

s  G  p{Ag).  This,  in  particular,  means  that  Tgg{s)  is  a 
compact  operator  and  thus  its  eigenvalues  A®  (s)  accumulate 
at  the  origin  as  |n|  — >  00  [13].  As  a  matter  of  fact  one  can 
make  a  much  stronger  statement. 

Lemma  4:  The  eigenvalues  A®  (s)  converge  to  the  origin 
uniformly  on 

Proof:  See  Appendix.  ■ 

Take  the  positive  integer  to  be  such  that  |A®  (s)|  <  e, 
s  G  D®,  for  all  |  n\  >  Nf.  Let  us  rewrite  (5)  as 

Z  det[X  -I-  ^Tgg{s)]  = 

^  n  (1  +  7A^(s))  +  Z  (H-7A^(s))  = 

|n|<tVe  |n|>Af5 

^  Z  (1 -f  7A®  (s))  +  ^  Z  (1 -I- 7A®  (s)).  (6) 

|n|<JV,  |n|>JV, 

It  is  clear  that  if  I7I  <  f  then  for  \n\  >  iV„  |7A«(s)|  < 
1,  and  1  +  7A®  (s)  can  never  circle  the  origin  as  s  travels 
around  S)®.  Thus  for  I7I  <  ^  the  final  sum  in  (6)  will  not 
contribute  to  the  encirclements  of  the  origin,  and  hence  we 


lose  nothing  by  considering  only  the  first  eigenvalues. 
There  still  remain  some  minor  technicalities. 

First,  let  denote  the  disk  |s|  <  e  in  the  complex 
plane.  Then  said  truncation  may  result  in  some  eigenloci 
(part  of  which  resides  inside  D^)  not  forming  closed  loops. 
But  notice  that  these  can  be  arbitrarily  closed  inside  D^,  as 
this  does  not  affect  the  encirclements  [2]. 

The  second  issue  is  that  for  some  values  of  s  G 
!Fgg{s)  may  have  multiple  eigenvalues,  and  hence  there 
is  ambiguity  in  how  the  eigenloci  of  the  Nyquist  diagram 
should  be  indexed.  But  this  poses  no  problem  as  far  as 
counting  the  encirclements  is  concerned,  and  it  is  always 
possible  to  find  such  an  indexing;  for  a  detailed  treatment 
see  [1]. 

Let  us  denote  by  A®  the  indexed  eigenloci  that  make  up 
the  generalized  Nyquist  diagram.  From  (6)  and  the  above 
discussion  it  follows  that 

G(o;det[X  +  7.F0,(s)]|^^^«)  = 

|n|<7V,  ^ 

which  together  with  Theorem  3  gives 

Theorem  5:  Assume  denotes  the  number  of  eigenval¬ 
ues  of  Ag  inside  C'*'.  For  and  as  defined  previously, 
Sp(Zl^^‘)  C  C-  for  I7I  <  i  iff 

~~  i  (An)|„|<7V7 

and 

(b)  E|„|<tv,  C'(- ->A^(s))  =pt- 

Finally,  the  closed-loop  system  G‘^*  is  exponentially  stable  iff 
Sp(Zl^')  C  C-  for  all  9  G  [0,  Ll].  ■ 

C.  Finite  Truncations  of  System  Operators 

The  above  development  means  that  for  a  given  e  >  0,  the 
eigenloci  that  fall  within  the  disk  =  {s  s.t.  |s|  <  e}  play 
no  role  in  the  Nyquist  stability  analysis  and  can  be  ignored 
as  long  as  |A|  >  e. 

This  suggests  that  one  could  truncate  Tgg,  or  equaiva- 
lently  truncate  Ag,  Bg,  Cg,  and  T,  and  effectively  treat  the 
stability  problem  as  one  for  a  family  of  (finite-dimentional) 
multivariable  systems  parameterized  by  the  variable  9. 

The  complication  here  is  that  although  a  truncation  re¬ 
moves  the  infinite  number  of  eigenloci  that  shrink  to  zero, 
it  also  affects  all  other  eigenloci,  no  matter  how  large  the 
truncation  is  taken  to  be.  Nevertheless,  it  can  be  shown  that 
by  increasing  the  size  of  the  truncation,  the  eigenloci  of  Tgg 
can  be  recovered  to  any  accuracy. 

Let  IlAr  be  the  projection  on  the  first  2N + 1  standard  basis 
elements  of  {c-n,  •  •  •  ,  cq,  •  •  •  ,  bat}.  Thus  IIjviFfJenjv  is 
the  {2N  -f  1)  X  {2N  +  1)  truncation  of  Tgg.  We  have 
Lemma  6:  If  C  G  piAN^FggA-N)  and 

WTgg-ANTggANW  ||(c-n^^g,n^)-i||  <  1, 

then  C  G  p{Tgg). 

Proof:  This  is  a  direct  consequence  of  Theorem  3.17, 
Chap  IV  of  [12].  ■ 

The  decay  of  the  diagonal  elements  of  gg,  and  the  decay 
with  increasing  I  of  the  //  elements  of  T,  can  be  used  to 


show  that  \\!FQ0—IIn!FG9^n\\  can  be  made  arbitrarily  small 
for  large  enough  N?  Hence,  by  Lemma  6  the  eigenloci  of 
^N^Qe^N  approximate  those  of  TQg  arbitrarily  well  as  N 
grows. 

D.  Regularity  in  the  9  Parameter 

Regarding  Theorems  3  and  5,  one  would  hope  for  some 
kind  of  ‘regularity’  with  respect  to  the  variable  9.  More 
precisely,  in  practice  one  would  like  to  plot  the  eigenloci 
and  check  the  Nyquist  stability  criterion  for  a  finite  number 
of  9i,  say  0i,  •  •  •  ,9^,  and  be  able  to  conclude  stability  for 
all  9  G  [0,  H]  if  the  9i  are  chosen  close  enough  to  each  other. 

It  is  possible  to  show  that  under  certain  mild  conditions 
on  A,  B  and  C,  all  points  of  the  plot  dei\I+"iTGo{s)\\^^^g 
change  continuously  with  9.  Moreover  we  can  show  that  the 
eigenloci  of  !FQe{s)  change  continuously  with  9.  We  do 
not  present  the  details  here. 

VI.  An  Illustrative  Example 

In  this  section,  we  will  consider  the  example  of  an  open- 
loop  system  governed  by 

dttp{t,x)  =  dl'ip{t,x)  +  i’{t,x)  +  u{t,x), 

y{t,x)  =  'tp{t,x), 

put  in  feedback  with  yF{x)  =  7C0s(a;).  Here  H  =  1,  and  B, 
C  are  the  identity  operator.  The  representation  of  the  system 
in  the  frequency  domain  is  dt'tpg{t)  =  (Ag  —  'j!F)ipg{t), 
where  Ag  =  diag{-  •  •  ,  —{9  +  n)^  +  1,  •  •  •  }  for  every  9  G 

[0, 1],  and  T  has  the  form  shown  at  the  end  of  Section  II 

with  fi  =  f-i  =  1  and  fi  =  0,  i  ^  ±1.  Notice  that  the 
open-loop  system  is  unstable. 

Recall  that  to  test  the  stability  of  the  closed-loop  system, 
one  has  to  apply  the  Nyquist  criterion  for  every  9  G  [0, 1]. 
We  take  Nyquist  paths  2)^  of  the  form  shown  in  Figure  4(a) 
with  r  =  20.  The  indentation  is  chosen  appropriately  as  to 
avoid  the  eigenvalues  of  Ag  at  the  origin.  Let  us  take  a  look 
at  the  Nyquist  plots  for  two  particular  values  of  9\ 

6*  =  0  :  A  =  0,0,1  are  the  eigenvalues  of  Ao  inside 
!D°,  hence  =  3,  and  we  need  three  counter-clockwise 
encirclements  of  — 1/7  to  achieve  closed-loop  stability.  As 
can  be  seen  in  Figure  4(b)  and  its  blown-up  version  (c), 
one  possible  choice  would  be  to  take  — 1/7  to  be  purely 
imaginary  and  — 0.2j  <  — 1/7  <  0.2j.  Clearly  such  — 1/7 
is  encircled  three  times  by  the  eigenloci. 

9  =  0.5  :  A  =  0.75,0.75  are  the  eigenvalues  of  AI0.5 
inside  hence  =  2,  and  we  need  two  counter¬ 

clockwise  encirclements  of  — 1/7  to  achieve  closed-loop 
stability.  Again,  from  Figure  4(d),  if  — 1/7  is  taken  to  be 
purely  imaginary  and  — 0.2j  <  — 1/7  <  0.2j,  then  — 1/7  is 
encircled  twice  by  the  eigenloci. 

As  a  matter  of  fact,  it  can  be  shown  that  — 0.2j  <  -1/7  < 
0.2j  stabilizes  the  closed-loop  system  for  all  values  of  9  G 
[0,1]. 

^This  is  also  a  direct  consequence  of  the  fact  that  TQg  G  ^i(f^)  is  a 
compact  operator  and  can  therefore  be  approximated  arbitrarily  well  by  a 
sequence  of  finite-dimensional  operators  [13]. 


Eig.  4.  (a)  The  Nyquist  path  X>®;  (b)  Nyquist  plot  for  6  =  0;  (c)  Blown-up 
version  of  plot  in  (b);  (d)  Nyquist  plot  for  d  =  0.5. 


VII.  Conclusions  and  Future  Work 

We  develop  a  generalized  Nyquist  stability  criterion  that 
is  applicable  to  a  class  of  spatially-distributed  systems  with 
infinite-dimensional  input/output  spaces  and  unbounded  sys¬ 
tem  operators. 

Future  work  in  this  direction  would  include  considering 
a  wider  class  of  spatially-distributed  systems,  for  example 
those  for  which  the  Fourier  symbol  of  A  is  itself  a  matrix 
(i.e.,  the  state  has  Euclidean  dimension  larger  than  one). 
Also,  the  affect  on  stability  of  different  basic  frequencies  H 
of  the  periodic  feedback,  and  the  occurrence  of  parametric 


resonance,  can  be  investigated  in  this  framework.  It  would  be 
of  interest  to  compare  such  results  to  those  already  obtained 
using  norm  and  spectral  analysis  [4]  [14]. 

VIII.  Appendix 

Proof  of  Lemma  1:  Since  s  G  p{Ag),  then  s  ^  A{6  +  Lln) 
for  any  n  G  Z,  which  together  with  (3)  and  the  boundedness 
of  B  and  C  implies  that  all  diagonal  elements  of  fle(s)  are 
finite.  Let  n  =  1,  2,  •  •  • ,  be  a  reordering  of  these 

elements  such  that  |pi(fJe(s))|  >  |p2(0e(s))|  >  •  •  •  >  0. 
Then  clearly  cr„(be(s))  =  |p„  (f/0(s))  |.  But  from 

Assumption  (i)  it  follows  that  ‘^n(Se(s))  <  oo.  Hence 

e^e(s)  G  and  since  T  G  TQe^s)  G 

[10].  ■ 


To  prove  Theorem  2  we  need  the  following  lemma. 

Lemma  Al:  For  s  G  p{Ag),  det[X  + 7iF^g(s)]  is  analytic 
in  both  7  and  s. 

Proof:  For  s  G  p{Ae),  'jTQe{s)  G  by  Lemma 

1.  Also  'jTQe{s)  =  ^TC0{sI  —  A0)~^Be  is  clearly  analytic 
in  both  7  and  s  for  s  G  p{Ag).  Then  it  follows  from  [10, 
pl63]  that  det[X+ 7iF5e(s)]  too  is  analytic  in  both  7  and  s 
for  s  G  p{Ag).  ■ 


Proof  of  Theorem  2:  Since  S®  does  not  pass  through 
any  eigenvalues  of  Ae,  s  G  p{Aq)  and  thus  ')TCq{sX  — 
Ae)~^Be  G  from  Lemma  1.  Then  from  [9],  (X  + 

')!FCe{sI  —  Ag)~^Bg)  ^  exists  and  belongs  to  iff 

det[X  +  ^TCe{sJ  —  Ae)~^Bg]  f  0,  which  is  satisfied 
by  assumption.  Applying  the  matrix  inversion  lemma  to 
{X  +  ^!FCg{sX  —  Ae)~^ Be)  we  conclude  that  s  G 

and  (sX  -  Af)-^  =  (sX  -  Ag  +  Bg-fTCg)-^  G 

Now  since  (sX  —  Ag)~^,  Bg,  Cg,  and  T  all  belong 
to  \  and  (sX  —  Ag)~^  G  we  have  (sX  — 

Ag)~^BgjBCe{sX  —  Ag)~^  G  [10].  Thus,  from  the 

identity 

{sX  -  Af)-^  -  {sX  -  Ag)-^  = 

-  (sX  -  Atr^BgjBCgisX  -  Ag)-\ 


it  follows  that  (sX  —  A‘g)~^  G  In  particular  (sX  — 

Ag)~^  and  (sX  —  Af)~^  are  both  in  3§o{£'^),  which  means 
that  Ag  and  Af  both  have  discrete  spectrum  [12]  and  in 

{sX  -  Afy^ds  -  ^  [  {sX-Ag)~^ds  = 

27rjJ2,e 

-^f  (sX-  At)-^BgjBCe{sX  -  Ag)-^ds 

each  term  on  the  left  side  is  a  finite-dimensional  projection 
[10,  pll,  pi 5].  Taking  the  trace  of  both  sides  and  changing 
the  order  of  integration  and  trace  on  the  right^  we  have 


tr  /  (sX  —  Af)  ^ds  —  tr  f  {sX  —  Ag)  ^ds 

.27rj7j)»  J  L27rj72je 

tv[{sX  -  At)-^B0-/BCe{sX  -  Ae)-^]ds.  (Al) 

jJsiO 


1 


^Since  (sX  — ^Bo'yJ^Ce(sT  —  Ae)  ^  G  =^i (•^^),  its  trace  is  well- 
defined  and  finite. 


From  [12],  the  left  side  of  (Al)  is  equal  to  the  number  of 
eigenvalues  of  As  in  c+  minus  the  number  of  eigenvalues 
of  Ag  in  C+. 

On  the  other  hand,  let  the  path  C®  be  that  traversed  by 
det[X  +  "f!FQg{s)\  as  s  travels  once  around  D®  (where  2)® 
lies  entirely  in  p{Ag)).  By  Lemma  Al,  dei\X +  ^!FQg{s)\  is 
analytic  in  s,  and  if  det[X  +  ')!FQg{s)\  7^  0  on  2)^  we  have 


C(0;det|l  +  ,:F5.W]|_^^,)  = 

i^det[X  +  7X-gg(5)]^^^ 
2'Kj  J^e  det[X  + 7X'(le(s)] 


(A2) 


Using  [10,  pl63]  and  the  assumption  dei\X+^!FQg{s)\  A 
s  G  23®,  we  arrive  at 


Adet[X  +  7X'gg(s)] 
det[X  +  '-fTQe{s)\ 


tr 


{X  +  ^Tgg{s)) 


=  -tr[(sX  -  Af)-^BgjBCg{sX  -  Ag)-^]. 


This,  together  with  (A2)  and  (Al)  gives  the  required  result.  ■ 


Proof  of  Lemma  4\  For  s  G  p{Ag),  dei\X  +  ^TQg{s)\ 
is  analytic  in  both  7  and  s  by  Lemma  AL  The  proof  now 
proceeds  exactly  as  in  [2,  pl40]  and  is  omitted.  ■ 
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